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Most animal-microbe symbiotic interactions must be advantageous to the host and provide nutritional benefits 
to the endosymbiont. When the host provides nutrients, it can gain the capacity to control the interaction, pro¬ 
mote self-growth, and increase its fitness. Chlorella-like green algae engage in symbiotic relationships with cer¬ 
tain protozoans, a partnership that significantly impacts the physiology of both organisms. Consequently, it is 
often challenging to grow axenic Chlorella cultures after isolation from the host because they are nutrient fastid¬ 
ious and often susceptible to virus infection. We hypothesize that the establishment of a symbiotic relationship 
resulted in natural selection for nutritional and metabolic traits that differentiate symbiotic algae from their 
free-living counterparts. Here, we compare metabolic capabilities of 5 symbiotic and 4 free-living Chlorella strains 
by determining growth levels on combinations of nitrogen and carbon sources. Data analysis by hierarchical clus¬ 
tering revealed clear separation of the symbiotic and free-living Chlorella into two distinct clades. Symbiotic algae 
did not grow on nitrate but did grow on two symbiont-specific amino acids (Asn and Ser) on which the free- 
living strains did not grow. The use of these amino acids was exclusively affected by the presence/absence of 
Ca 2+ in the medium, and the differences were magnified if galactose was provided rather than sucrose or glu¬ 
cose. In addition, Chlorella variabilis NC64A genomic and differential expression analysis confirmed the presence 
of abundant amino acid transporter protein motifs, some of which were expressed constitutively both axenically 
and within the host. Significantly, all 5 symbiotic strains exhibited similar physiological phenotypes even though 
they were isolated as symbionts from different host organisms. Such similarities indicate a parallel coevolution of 
shared metabolic pathways across multiple independent symbiotic events. Collectively, our results suggest that 
physiological changes drive the Chlorella symbiotic phenotype and contribute to their natural fitness. 

© 2016 Elsevier B.V. All rights reserved. 


1. Introduction 

Successful endosymbiosis provides advantages to both the host and 
the endosymbiont. Benefits may include better adaptation to limited 
nutrients or reduction of mortality via protection against damage by 
UV light or pathogens (e.g., viruses). In such scenarios, symbionts in¬ 
crease their reproductive capacity and fitness within their hosts relative 
to possible non-host environments [3,14,18,19]. For example, some 
protozoans harbor intracellular green algae in an inherited mutually 
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beneficial symbiotic relationship, which serves as a well-recognized 
model for studying endosymbiotic relationships [27,28,37,45]. These 
unicellular Chlorella -like green algae, often referred to as zoochlorellae, 
inhabit the gastrodermal symbiosomes (perialgal vacuoles) of different 
protozoans, and transfer a significant amount of their photosynthetical- 
ly fixed carbon (e.g., maltose, fructose) to their non-photosynthetic 
partners [7,18,31]. In this context, symbiotic Chlorella spp. still require 
nutrients such as nitrogen, which they obtain from the host and then as¬ 
similate into the algal metabolome [33,50]. The mechanisms involved in 
these interactions have not been completely elucidated; however, the 
metabolic pathways involved in nitrogen (N) and carbon (C) utilization 
could be crucial physiological signatures of the endosymbiotic existence 
[32]. Therefore, elucidating how such processes work would open new 
avenues of research in the understanding of the molecular, cellular, and 
organismal adaptations that allow successful mutualism. 


http ://dx.doi.org/l 0.1016/j .algal.2016.06.001 
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Protozoa -Chlorella interactions can be disrupted, and some attempts 
to isolate intact algae free of the host have been successful. These in¬ 
clude algae that associate with several species of protozoans, including 
Paramecium bursaria [23,45], Acanthocystis turfacea [23], and Hydra 
viridis [32]. Another approach for identifying ex-symbiotic algal strains 
has relied on their susceptibility to large DNA virus infections after dis¬ 
ruption of the host-Chlorella interaction [22,30,34,47]. The only docu¬ 
mented symbiotic, virus-susceptible Chlorella strains that have been 
cultured axenically include Chlorella variabilis NC64A 1 [49], C. variabilis 
Syngen2-3 [47], C. variabilis F36-ZK [11,16,38], C. variabilis OKl-ZK ([11, 
38], Quispe et al., manuscript in preparation), and C. heliozoae SAG 3.83 
[6]. For the purpose of this paper, these symbiotic virus-susceptible algal 
strains will be referred to as symbiotic algae. 

We have studied the Chlorella-v irus interaction for the past 35 years 
and have been aware of the fastidious nutrient requirements possessed 
by symbiotic algae [1,20,33]. For example, unlike most Chlorella species, 
the symbiotic algae do not grow on Bolds' Basal Medium (BBM), which 
has nitrate (N0 3 ) as its sole N source. Consequently, 0.1% peptone is 
added to BBM for axenic growth of these symbiotic strains [15,16,25, 
47,49]. We hypothesize that this requirement reflects a past symbiotic 
relationship that spurred selection for specific nutritional and metabolic 
features present in symbiotic algae. In this study, we test this idea by an¬ 
alyzing some physiological traits and growth requirements in 4 free- 
living and 5 symbiotic Chlorella species. 

Our physiological evaluation focused on alternative N and C sources. 
The results indicate that symbiotic algae are better able to assimilate N 
and C sources not normally available to the free-living strains. Signifi¬ 
cantly, they prefer organic N sources rather than the inorganic N sources 
(e.g., N0 3 or NH 4 ), which are the primary N sources in the environment. 
Importantly, all symbiotic strains tested exhibit similar metabolic phe¬ 
notypes even though they are polyphyletic and may have arisen as pro¬ 
tozoan symbionts from several independent symbiotic events [13,38]. 
Importantly, these similarities denote a parallel coevolution of similar 
metabolic pathways across multiple independent symbiotic events. 
Taken together, this evolutionary genome plasticity and metabolic reg¬ 
ulatory rewiring could come at a cost in the form of the inability of sym¬ 
biotic Chlorella to survive as free-living organisms in virus replete and/ 
or nutrient limiting environments [39]. 

2. Materials and methods 

2.1. Algal strains 

Symbiotic C. variabilis ATCC 50258 (NC64A), C. variabilis ATCC 30562 
(Syngen 2-3), and C. heliozoae SAG 3.83 (SAG 3.83) were maintained as 
slant stocks at 4 °C. Symbiotic C. variabilis NIES-2540 (F36-ZK) and 
C. variabilis NIES-2541 (OKl-ZK) were obtained from the Japanese Cul¬ 
ture Collection of the National Institute for Environmental Studies 
(http://mcc.nies.go.jp). Stock samples of free-living Chlorella strains 
Chlorella sorokoniana UTEX-1230 (UTEX-1230), Cyamus kessleri UTEX- 
2228 (B228), and Chlorella protothecoides UTEX-29 (CP-29) were ob¬ 
tained from the Culture Collection of Algae at the University of Texas 
at Austin (https://utex.org), and C. sorokoniana CCTCC M209220 (CS- 
01) was obtained from Minxi Wan at Johns Hopkins University. The se¬ 
lection for free-living strains was based on the proposed phylogeny of 
Chlorella species published by Rosenberg et al. [41]), from which we 
chose four representative strains. 

2.2. Cell cultures 

Symbiotic and free-living strains were grown on BBM [3] supple¬ 
mented with 0.1% (w/v) peptone, 0.5% (w/v) sucrose, and 0.001% (w/v) 
thiamine (complete MBBM). All of the BBM modifications had N0 3 and 
sucrose omitted (N-/C-BBM) and being replaced by the labeled N and C 
sources. Where indicated, 0.1% peptone was replaced with 0.1% (w/v) 
casamino acids. The ability of algae to use different N and C sources 


was tested by adding them to unsupplemented BBM (N-/C-BBM). 
Thus, 0.22 pm filter-sterilized stock solutions of N and C sources were 
added to a final concentration of 10 mM. All flasks were supplemented 
with 0.001% (w/v) thiamine. To test the effect of Ca 2+ deprivation on 
algal growth, we used a C-, N-, and Ca 2+ -deficient BBM (N-/C-/Ca 2+ - 
BBM), i.e., Ca 2+ was not included in the BBM. 

The algae were grown in 125 ml narrow mouth Erlenmeyer flasks 
with 30 ml of supplemented BBM. For the inoculum, MBBM log-phase 
actively growing cells were pelleted and washed 3 times with either 
N-/C-BBM or N-/C-/Ca 2+ -BBM medium. Flasks were inoculated to a 
final low cell density of 3-5 x 10 4 cells/ml and shaken at 26 °C and 
180 rpm in continuous light for variable time periods because the sym¬ 
biotic growth rates were slower than their free-living counterparts. 
Free-living strains were grown for 9 days on BBM with added N or for 
7 days when both N and C were added. Similarly, symbiotic strains 
were grown for 12 days on BBM with added N or for 9 days when 
both N and C were included. MBBM and un-supplemented BBM were 
used as controls. Triplicate samples were used for the symbiotic algae, 
and duplicate samples were used for the free-living strains. Photographs 
of the flasks were taken with a 12.1 M pixel Sony Cyber-shot digital 
camera and organized using Adobe Photoshop CS5.1. They are shown 
in Supplementary Figs. 1-6. 

2.3. Hierarchical clustering analysis 

Cluster 3.0 for Mac OS X (http://rana.lbl.gov/EisenSoftware.htm) 
and JavaTreeView Version 1.1.6r4 (http://jtreeview.sourceforge.net/) 
programs were used to analyze and quantify the growth experiments. 
The hierarchical clustering algorithm was performed using the 
average-linkage method applied to the data set. This algorithm pro¬ 
duced a dendrogram that assembled all elements into a single tree, 
which arranged the strains and treatments according to similarities in 
their growth patterns. The data set consisted of rows and columns 
representing the 9 algal strains and the numerical score for each 
media condition. Analyses were performed both on the bulk data and 
as subsets by treatment. Numerical scores were assessed for individual 
flasks using a 0 to 5 scale, with 5 representing the best growth and 0 
the absence of visible growth. The data sets were represented graphical¬ 
ly in hierarchical clusters by coloring each cell on the basis of the numer¬ 
ical flask score. Flasks with scores of 0 were colored black while the 
higher scores were reds of increasing intensity to denote growth. The 
dendrogram is attached on both axes to the colored graph to indicate 
the computed relationships among both growth conditions and Chlorel¬ 
la species. 

2.4. Comparative genomics 

Members of a collection of characterized AA transporters from 
Arabidopsis thaliana [46] were used to perform reciprocal BLAST 
searches [2] against C. variabilis NC64A and C. sorokoniana UTEX-1230 
(UNL algal consortium, in preparation) genomes, using a value of 
1 x 10“ 10 as a cutoff. Each algal protein identified an A. thaliana AA 
transporter, and the gene designations and E-values for each gene are 
presented in Supplementary Tables 1 and 2, respectively. Similarly, 35 
putative AA transporters from NC64A [4] were used to perform a recip¬ 
rocal BLAST search against the UTEX 1230 proteome using a value of 
1 x 10“ 10 as a cutoff (Table 1). The orthology of candidate sequences 
was verified using the KEGG database [17]. 

2.5. RNAseq analysis 

Data sets from RNAseq experiments were downloaded to the public 
Galaxy platform server (www.usegalaxy.org) and manipulated with 
data analysis tools as described below. Tophat settings were as defined 
by the usegalaxy.org defaults (Galaxy Version 0.9) except that the max¬ 
imum intron length was set at 1500 bp. For axenic C. variabilis NC64A 
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Table 1 

Accession and scaffold numbers of putative C. variabilis NC64A orthologs to Chlorella sorokiniana UTEX-1230 proteins involved in AA transport. 



NC64A gene ID 

Protein ID 

Scaffold in C. sorokoniana 

e-value 

Bit score 

Identity 

Protein length (aa) 

C. variabilis C. sorokoniana 

1 

138810 

EFN52376 

18.g93.isol 

3E-073 

231 bits (590) 

126/181 (70%) 

183 

504 

2 

50436 

EFN58616 

4.gl63.isol 

5E-088 

284 bits (726) 

186/451 (41%) 

410 

690 

3 

138809 

EFN52375 

18.g93.isol 

2E-113 

340 bits (871) 

190/306 (62%) 

287 

504 

4 

144770 b 

EFN56324 

56.g4.iso2 

7E-128 

380 bits (977) 

193/361 (53%) 

471 

405 

5 

144819 

EFN56345 

4.g37.iso3 

2E-120 

364 bits (935) 

213/389 (55%) 

489 

469 

6 

142340 a 

EFN58068 

34.g94.isol 

IE-134 

426 bits (1096) 

219/305 (72%) 

695 

1070 

7 

37093 ab 

EFN51991 

172.gl06.isol 

2E-141 

431 bits (1108) 

244/446 (55%) 

519 

815 

8 

135113 a 

EFN54610 

170.g37.isol 

3E-134 

414 bits (1063) 

246/309 (80%) 

932 

461 

9 

145403 

EFN55845 

4.gl63.isol 

2E-134 

410 bits (1053) 

255/547 (47%) 

535 

690 

10 

133029 

EFN59609 

28.gl35.isol 

5E-139 

409 bits (1052) 

257/443 (58%) 

431 

449 

11 

134730 

EFN54961 

57.g52.isol 

2E-156 

455 bits (1171) 

261/424 (62%) 

453 

469 

12 

134234 a 

EFN55146 

33.gl85.isol 

4E-142 

420 bits (1079) 

263/484 (54%) 

473 

480 

13 

142091 

EFN57962 

8.gll6.iso2 

2E-158 

461 bits (1186) 

266/417 (64%) 

471 

466 

14 

135437 3 

EFN54731 

3.g27.isol 

4E-161 

475 bits (1222) 

272/399 (68%) 

518 

606 

15 

51413 

EFN57306 

84.gl38.isol 

IE-174 

503 bits (1294) 

282/452 (62%) 

452 

490 

16 

49669 

EFN60144 

4.gl58.isol 

8E-179 

533 bits (1372) 

294/470 (63%) 

726 

742 

17 

133360 a 

EFN60071 

174.g65.isol 

0 

552 bits (1422) 

296/405 (73%) 

692 

685 

18 

56488 a 

EFN59984 

53.g43.isol 

0 

585 bits (1509) 

300/495 (61%) 

973 

515 

19 

138717 

EFN52627 

21.g78.isol 

0 

547 bits (1410) 

317/540(59%) 

498 

547 

20 

133351 

EFN60067 

27.g54.isol 

0 

600 bits (1548) 

331/494 (67%) 

489 

473 

21 

138133 a,b 

EFN59501 

110.g43.isol 

0 

594 bits (1532) 

336/587 (57%) 

576 

512 

22 

58128 a,b 

EFN54604 

15.gl50.isol 

0 

637 bits (1644) 

337/516(65%) 

522 

484 

23 

57473 a ' b 

EFN56726 

91.g67.iso3 

0 

681 bits (1758) 

346/471 (73%) 

476 

468 

24 

138505 

EFN52920 

43.g21.isol 

0 

659 bits (1701) 

353/546 (65% 

544 

560 

25 

32765 a,b,c 

EFN51990 

35.gll7.iso2 

0 

659 bits (1701) 

366/650 (56%) 

605 

635 

26 

59057 a ' b 

EFN51898 

35.gll7.iso2 

0 

645 bits (1664) 

373/698 (53%) 

742 

635 

27 

133392 

EFN60084 

39.g82.isol 

0 

688 bits (1775) 

378/516 (73%) 

516 

496 

28 

59479 a,b,c 

EFN50622 

172.gl06.isol 

3E-056 

194 bits (492) 

97/158 (61%) 

270 

815 

29 

133952 

EFN55688 

28.gl05.isol 

2E-023 

101 bits (251) 

50/88 (57%) 

336 

674 

30 

133953 

EFN55689 

28.gl02.isol 

7E-053 

172 bits (435) 

93/183 (51%) 

197 

238 

31 

24724 a,b 

EFN54340 

76.g4.isol 

IE-051 

177 bits (450) 

107/269 (40%) 

227 

560 

32 

53357 b 

EFN53996 

34.gl91.isol 

4E-066 

224 bits (572) 

166/486 (34%) 

385 

689 

33 

137496 a 

EFN53131 

27.g54.isol 

7E-047 

176 bits (446) 

144/415 (35%) 

1010 

473 

34 

138482 

EFN52813 

177.gll.isol 

IE-065 

220 bits (560) 

215/445 (48%) 

425 

489 

35 

17797 b 

EFN50713 

172.g67.isol 

4E-031 

115 bits (287) 

54/87 (62%) 

92 

815 

36 

142334 b 

EFN58316.1 

13.g249.iso3 

0 

661 bits (1706) 

365/420 (87%) 

535 

487 

37 

140447 b 

EFN58455.1 

106.g215.iso2 

1.00E-96 

302 bits (774) 

142/179 (79%) 

376 

576 

38 

7483 b 

EFN50706.1 

92.gll6.iso2 

0 

667 bits (1720) 

373/497 (75%) 

493 

624 

39 

58448 b 

EFN53780.1 

136.g46.isol 

0 

560 bits (1442) 

316/671 (47%) 

639 

1388 

40 

36103 b 

EFN54400.1 

139.gl41.isol 

2E-130 

400 bits (1028) 

247/427 (58%) 

444 

832 


a C. variabilis gene ID of expressed AA transporter genes during axenic growth ([4] and [42]), see Fig. 7. 

b C. variabilis gene ID of differentially expressed AA transporter genes during axenic and endosymbiont states ([42] and [26]), see Table 2. 
c Genes (32765 and 59479) that accounted for the majority of the mapped differentially expressed reads in axenic and symbiont growth. 


(Fig. 7), we used an uninfected control sample (NCBI SRA accession 
SRX316780) from a published viral infection experiment conducted in 
our lab [42]. For C. variabilis growing endosymbiotically within 
P. bursaria, we downloaded RNAseq data sets (NCBI SRA accessions 
DRX003053, DRX003054, and DRX003055) [26]. These sequence files 
were reported to contain RNAseq reads mapping to the NC64A genome, 
thus providing potential information regarding algal genes that are dif¬ 
ferentially expressed when grown axenically on MBBM versus its natu¬ 
ral endosymbiont stage within P. bursaria. 

The FASTQ files were converted to FASTQSANGER format with the 
FASTQGroomer tool (default settings) [5] and then Tophat [24] was 
used to align these data sets to the NC64A genome assembly [4] with 
minimum and maximum intron lengths of 50 and 5000, respectively. 
Around 1% (-970,000) of the P. bursaria derived reads aligned to the 
C. variabilis NC64A genome, and these reads were taken to represent 
a snapshot of gene expression in the endosymbiont cells. The same 
analysis pipeline was applied to the axenic C. variabilis NC64A data. 
Reads that mapped to the genomic intervals for each putative AA 
transporter in Supplementary Table 1 were counted using the Inte¬ 
grated Genome Browser software package [36] and normalized as 
total mapped reads per gene in each condition per million mapped 
reads. 


3. Results 

3 A. An overview of nitrogen metabolism in symbiotic algae 

This study addresses difficulties in the isolation and growth of sym¬ 
biotic algae. For example, C. variabilis NC64A grows well on BBM with 
0.1% peptone and 0.5% sucrose (MBBM) but does not grow on un¬ 
supplemented BBM [15,16,25,47,49]. Additionally, previous reports 
established that peptone supplemented BBM was sufficient for growth 
of C. variabilis NC64A and C. variabilis F36-ZK [16]. We have now includ¬ 
ed three more symbiotic strains (C. variabilis OKl-ZK and Syngen 2-3 as 
well as C. heliozoae SAG 3.83), and these 5 Chlorella strains comprise the 
symbiotic algae used in this study. Our focus on N metabolism in the 
symbiotic algae was prompted by our observation that NC64A grew al¬ 
most as well in a defined galactose-urea medium (BBM with 10 mM ga¬ 
lactose and 10 mM urea) as it did in MBBM. Two of the other symbiotic 
algae (SAG 3.83 and Syngen 2-3) also grew well when N0 3 and sucrose 
in BBM were replaced with 10 mM galactose and 10 mM urea, while 
F36-ZK and OKl-ZK grew well with the addition of 0.001% thiamine. 
Since galactose-urea-BBM is a chemically defined growth medium, we 
conclude that three of the symbiotic algae have no vitamin require¬ 
ments while two require thiamine. In particular, we did not find 
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evidence that symbiotic Chlorella required biotin or cobalamin/vitamin 
E>i 2 [9]. It is quite possible that these strains would grow better with 
added vitamin E>i 2 [23]. For convenience, 0.001% thiamine was added 
to all subsequent test media. 

3.2. Nutritional analyses identified physiologic signatures for symbiotic and 
free-living chlorella species 

The 5 symbiotic and 4 free-living Chlorella strains [(C. sorokoniana 
(UTEX-1230), C. sorokoniana (CS-01), C. kessleri (B228), and 
C. protothecoides (CP-29)] were compared with regard to their ability to 
utilize 64 different N and C sources. For each combination, a low cell den¬ 
sity of Chlorella (3-5 x 10 4 cells/ml) was inoculated into the indicated me¬ 
dium and cell growth was assessed by the intensity of the green color 
after 7-12 days (Supplementary Figs. 1-6). Triplicate and duplicate sam¬ 
ples were used for the symbiotic and free-living algae, respectively. The 
matrix of 64 growth conditions and 9 Chlorella strains was analyzed by 
two-way heat maps (Fig. 1). These average-linkage maps display differ¬ 
ences in the metabolic capabilities of the 5 symbiotic and 4 free-living 
Chlorella. The organic and inorganic N sources were tested with or with¬ 
out one of three sugars (sucrose, glucose, or galactose) as a C source. 
Both the C and N sources were routinely added at 10 mM. The rows in 
Fig. 1 represent combinations of 3C and 12 N sources (2 complex mix¬ 
tures, 7 organic, and 3 inorganic) or the N sources alone, while the col¬ 
umns represent the growth levels achieved by the 9 Chlorella strains on 
the 64 media combinations. The 9 Chlorella strains clearly separated into 
two clusters based on their nutritional capabilities: a symbiotic clade 
and a free-living clade (Fig. 1). Purple labels identify inorganic N sources 
present at 1 mM instead of 10 mM while orange labels represent media 
prepared without Ca 2+ . As a major difference, 4 out of the 5 symbiotic 
Chlorella did not grow on N0 3 as the sole N source, either with or without 
added sugar, whereas the 4 free-living Chlorella grew well on N0 3 with 
sugar and poorly on N0 3 without sugar (Fig. 1 ). 

3.3. Casamino acids and peptone supply all the C and N needed by symbiotic 
algae 

The heat map presented in Fig. 1 is cumbersome to read due to its 
large size. Therefore, we prepared 41 smaller analyses reflecting 


metabolic subsets of the data and of these the 6 most relevant subsets 
are presented here. In each subdivision, the 5 symbiotic and 4 free- 
living Chlorella were compared based on their growth on the respective 
supplemented BBMs. In the first case (Fig. 2), we analyzed peptone and 
casamino acids with and without sucrose. The symbiotic and free-living 
strains formed two distinct clades (Fig. 2) while two media clusters ap¬ 
peared, one for casamino acids alone and another for peptone alone, 
peptone and sucrose (MBBM), and casamino acids and sucrose. The 
free-living strains did not grow as robustly as the symbiotic Chlorella 
did on casamino acids alone (Fig. 2). 

All the symbiotic Chlorella grew better on 0.1% casamino acids than 
on 0.1% peptone and, in all cases removal of sucrose and N0 3 from the 
control MBBM had no effect on their growth (Fig. 2). Although the dif¬ 
ferences were slight, the results suggest that organic N sources rich in 
AAs (casamino acids) were better assimilated by the symbiotic Chlorella 
than the free-living strains. By contrast, the free-living Chlorella 
responded better to the addition of sucrose, suggesting that they were 
better adapted to the presence of a sugar source. 


3.4. Asn and Ser were better assimilated by symbiotic Chlorella 

The results in Fig. 2 suggested that free AAs were better assimilated 
by the symbiotic strains. Therefore, we tested the 9 Chlorella strains on 
10 organic and inorganic N sources. The dendrogram (Fig. 3) shows a 
clear separation between the symbiotic and free-living Chlorella based 
on their N assimilation patterns. All of the Chlorella grew with Arg, 
urea, Gin, and Gly as the sole N source but the symbiotic grew slightly 
better. Additionally, a separate cluster was formed by Asn and Ser 
(Fig. 3). Again, the symbiotic algae grew better on these AAs, with Asn 
being better than Ser (Fig. 3). Thus, Asn and Ser appeared to be 
symbiont-specific in that they were used poorly, if at all, by free-living 
Chlorella. Hence, Asn and Ser might be important N sources during sym¬ 
biotic growth. Pro also clustered with Asn and Ser (Fig. 3), but it did not 
support growth of the symbiotic Chlorella strains as well as Asn and Ser. 
It is important to note that these AAs were not simply fulfilling auxotro¬ 
phic requirements because all symbiotic Chlorella were fully prototro¬ 
phic for AA biosynthesis. That is, they all grew on minimal defined 
media with urea as the sole N source. 



Fig. 1 . Hierarchical heat map (average-linkage) clusters of symbiotic and free-living Chlorella strains based on their metabolic capabilities. Columns represent combinations of organic and 
inorganic nitrogen (N) sources with or without the addition of a carbon (C) source. All were added at 10 mM concentrations except purple labels identify inorganic N sources at 1 mM 
concentration, and orange labels represent media without Ca 2+ . MBBM medium was the control. Rows represent the 5 symbiotic strains (green) and 4 free-living strains (blue). Tree 
diagrams indicate the nature of the computed relationship among growth conditions and among Chlorella species. A color scale indicates relative growth: red represents robust growth 
and black represents absence of growth. Flask tests were performed for 9-12 days for the symbiotic and 7-9 days for the free-living strains. Subsequent heat map Figs. 2-6 follow a 
similar layout. 

















































336 


C.F. Quispe et al /Algal Research 18 (2016) 332-340 


2 

0 


CD 

C/5 


CO 

O 

k_ 


0 

o 

D 

CO 


CO 

+ 


o 

CO 

CO 

o 

"O 

"D 

o 

O 

3 

CO 

0 

0 

+ 

O 

O 

0 0 

c 

c 

c c 
o o 

-l-» -4— ' 

E 

0 

E 

0 

Q. Q. 

CO 

CO 

0 0 

0 

0 

Q- Q- O 

o 


Relative growth 



Max growth No growth 




SAG 3.83 
Syngen 2-3 
NC64A 
OK1-ZK 
F36-ZK 
CS-01 
UTEX 1230 
B228 
CP 29 


0 

0 

CD 

m 

E~ E 

3 C 3 

C E C 
0^0 
E:§ E CO 
E o E E 5 2 -E -2^ 5 q 

<C0<<3OO<C0 


2 

Q_ 


Relative growth 



n 



SAG 3.83 
Syngen 2-3 
NC64A 
F36-ZK 
OK1-ZK 
UTEX 1230 
CS-01 
CP 29 
B228 


Fig. 3. Heat map subgroup from Fig. 1 compares inorganic and organic N sources at 10 mM 
concentrations as the sole N source. Inorganic sources include NH 4 tartrate, sodium N0 3 , 
and NH 4 acetate. Organic sources include arginine (Arg), urea, glutamine (Gin), glycine 
(Gly), asparagine (Asn), serine (Ser), and proline (Pro). A color scale indicates relative 
growth. Flask tests were performed for 12 days for the symbiotic and 9 days for the 
free-living strains. 


Fig. 2. Heat map subgroup from Fig. 1 displays variations of MBBM (sucrose + peptone). 
Peptone is replaced by 1% casamino acids. A color scale indicates relative growth. Flask 
tests were performed for 9 days for the symbiotic and 7 days for the free-living strains. 


AA transport is coupled to movement of ions, including Na + , H + , K + , 
Ca 2+ , and/or Cl - as well as movement of sugars. Additionally, Kato and 
Imamura [21 ]) showed that the presence of divalent but not monova¬ 
lent cations decreased AA uptake in C. variabilis F36-ZK. We also noticed 
minimal growth of 3 symbiotic strains when Ca 2 + was present in the 
media. Therefore, the influence of Ca 2+ on AA transport and possible 
differences between symbiotic and free-living Chlorella species was 
also examined. We compared the effectiveness of organic N sources 
both with (Fig. 4, black color) and without Ca 2+ (Fig. 4, orange color). 
Ca 2 + is one of the salts normally present in BBM so the Ca 2 + -free medi¬ 
um was made by not adding Ca 2+ . The method ensures that trace levels 
of Ca 2+ will be presented from the water, inoculum, and other media 
components. The effect of Ca 2+ on the assimilation of organic N sources 
was variable. However, we observed that Ca 2+ influenced the assimila¬ 
tion of Asn, Ser, and Pro in symbiotic growth but it was strain specific. By 
contrast, the absence of Ca 2+ had no appreciable effects on Asn, Ser, and 
Pro uptake in the free-living Chlorella , but strain-specific differences 
were observed for urea, Arg, Gly, and Gin. We concluded that Ca 2 + 
ions influenced the assimilation of Asn, Ser, and Pro in symbiotic algae 
and the assimilation of Arg, Gly, and Gin in free-living Chlorella. 

3.5. Inorganic nitrogen sources 

The only inorganic N source that clustered within the organic group 
was NH 4 acetate (Fig. 3). NH 4 tartrate and sodium N0 3 clustered in a dif¬ 
ferent clade which, without added sugar, gave poor or no growth for all 
strains (Fig. 3). This finding was surprising because NH 4 and N0 3 are the 


primary sources of N in most aquatic and marine environments [43]. 
The symbiotic algae exhibited a remarkable duality in that they were 
unable to use extracellular NH 4 or N0 3 but were able to use the intracel¬ 
lular NH 4 formed after Arg, Gin, Asn, and urea uptake. These observa¬ 
tions confirm reports on C. variabilis F36-ZK regarding the loss of N0 3 
assimilation coupled with an enhanced ability to take up certain AAs 
[20]. F36-ZK uses a pH-dependent proton symport for general AA trans¬ 
port [16]. 
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Fig. 4. Heat map subgroup from Fig. 1 compares removal of Ca 2+ (orange) from media 
with organic N sources. Flask tests were performed for 12 days for the symbiotic and 
9 days for the free-living strains. 
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3.6. Low NH 4 concentrations are better assimilated by symbiotic strains 


Fungal growth media commonly include a 10 mM N source. Howev¬ 
er, 10 mM NH 4 could be inhibitory or toxic to some Chlorella strains 
(Fig. 3) and thus, we compared NH 4 and N0 3 at 1 and 10 mM in both 
the absence (Fig. 5) and in the presence of sugars (Fig. 6A). NH 4 tartrate 
was chosen because it does not acidify fungal growth media as much as 
(NH 4 ) 2 S0 4 or NH 4 C1 [29]. Lower levels of NH 4 acetate and NH 4 tartrate 
(1 mM) supported growth of all the symbiotic strains while higher con¬ 
centrations (10 mM) were inhibitory (Fig. 5). Similarly, Minaeva and 
Ermilova [35] demonstrated that in NH 4 -supplemented medium, 
C. variabilis NC64A showed significantly lower growth rates than in 
the peptone-containing MBBM. Thus, we conclude that the symbiotic 
Chlorella exhibit NH 4 toxicity. Both symbiotic and free-living strains 
grew better on NH 4 acetate (1 mM) than NH 4 tartrate (1 mM). This in¬ 
ability to utilize NH 4 efficiently might be different if we had included vi¬ 
tamin B 12 in our basal medium since Kessler and Huss [23]) reported 
that related Chlorella strains grow normally on NH 4 in a Bi 2 - 
containing medium. However, our current understanding of the meta¬ 
bolic function of Bi 2 in algae focuses on its role in methionine biosynthe¬ 
sis, not NH 4 utilization [9]. 

The lower N0 3 level (1 mM) did not support any algal growth while 3 
of the free-living and Syngen 2-3 (symbiotic group) had minimal growth 
at the higher level (10 mM) (Fig. 5). This is the only cluster analysis where 
a symbiotic strain did not group within the symbiotic cluster. A separate 
question concerns why the 4 free-living Chlorella only exhibit minimal 
growth on N0 3 (Fig. 5). The answer is likely a combination of inoculating 
at low cell density (3-5 x 10 4 cells/ml) and only following growth for the 
first 9 days. Additionally, all free-living Chlorella grew well on N0 3 when 
either glucose or sucrose was added (Fig. 6A). Presumably, the sugars pro¬ 
vided extra reducing power needed for the N0 3 to NH 4 conversion. In 
summary, symbiotic algae possess an efficient system to import and me¬ 
tabolize many AAs and small oligopeptides but they cannot efficiently uti¬ 
lize N0 3 or NH 4 as sole N sources. 
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Fig. 6. Heat map subgroups from Fig. 1. Inorganic (A) and organic (B) N sources 
supplemented with glucose, sucrose or galactose. Purple labels identify N sources at 
1 mM concentration. A color scale indicates relative growth. Flask tests were performed 
for 9 days for the symbiotic and 7 days for the free-living strains. 


3.7. Galactose is a better C source than glucose or sucrose for symbiotic 
Chlorella but not for free-living Chlorella 

All of the algal growth levels in Figs. 3 and 5 were lower than those 
observed on MBBM, which contains both sucrose and peptone. Thus, 
the presence or absence of sugars is important for an analysis of N me¬ 
tabolism. In this regard, Schlee and Komor [44]) showed that Chlorella 
had a single high-affinity NH 4 transporter whose synthesis was re¬ 
pressed by NH 4 but stimulated by glucose. Similarly, Cho and Komor 
[8]) showed that Chlorella vulgaris expressed multiple AA transport sys¬ 
tems, the synthesis of which was not repressed by NH 4 but most were 
stimulated by glucose. Accordingly, we supplemented all N sources 
with 3C sources: sucrose (present in MBBM), glucose, or galactose. Ga¬ 
lactose was chosen because in other microbes it does not exert catabo- 
lite repression [12]. Results for the inorganic and organic N sources are 
reported in Figs. 6A and 6B, respectively. They should be compared 
with Figs. 3 and 5, which use the same N sources but without added 
sugars. None of the symbiotic strains grew with any combination of 
sugar and N0 3 (Fig. 6A) except for Syngen 2-3. Thus, even after addition 
of sugars, 4 of the 5 symbiotic strains were unable to utilize N0 3 as a sole 
N source (Fig. 6A). Similarly, as had been observed in Fig. 5, 1 mM NH 4 
salts gave better growth than 10 mM regardless of the sugar used, and 
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better growth was achieved with 1 mM NH 4 acetate than with 1 mM 
NH 4 tartrate. 

In all cases, the symbiotic Chlorella grew with galactose or sucrose 
while they grew poorly with glucose (Figs. 6A and 6B). The symbiotic 
Chlorella grew as well in galactose and organic N source as they did in 
MBBM (Fig. 6B). NC64A appears to have more enzymes involved in car¬ 
bohydrate metabolism than other sequenced chlorophytes [4], includ¬ 
ing many which are related to galactose metabolism. Interestingly, 
added glucose inhibited the symbiotic Chlorella on most organic N 
sources, with the most glucose-sensitive strain being NC64A (Fig. 6B). 
By contrast, the free-living strains had similar or better growth with 
sugars than on MBBM, and they preferred sucrose and glucose 
(Fig. 6A). Galactose had inhibitory effects on the growth of most free- 
living strains with all of the organic N sources except Gly. Strain B228 
was the only free-living strain able to utilize multiple organic N sources 
in the presence of galactose (Fig. 6B). 

3.8. Bioinformatic and transcriptomic analysis of amino acid transporter 
orthologs in Chlorella species 

AA transporters act as extracellular and intracellular nutrient sen¬ 
sors as well as transporters in all domains of life [10]. The five symbionts 
studied must have efficient systems for importing AAs from the 
P. bursaria host, which is reflected in their continued ability to use 
some AAs as a source of N instead of inorganic N sources (Figs 2, 3, 
and 6B as well as [1,4,18-20]). Thus, the physiological observations 
depicted in Figs. 1-6 led us to hypothesize that in nature the protozoan 
host regulates the population of symbiotic Chlorella by using AAs to re¬ 
strict their N supply. Keeping the N supply low and the chlorophyll con¬ 
tent high (5- to 10-fold higher) is consistent with the symbiont 
functioning to provide excess photosynthate to the host in the form of 
secreted maltose [40]. It seemed likely that the function of symbiont- 
specific AA transporters might persist following their release from sym¬ 
biosis. In order to evaluate potential contributions of endosymbiosis on 
the evolution and expression of the AA transporters identified, we com¬ 
pared the genes for AA transporters in the symbiont NC64A and the 
free-living alga C. sorokoniana UTEX-1230 (Table 1 ), as well as analyzing 
gene expression profiles of NC64A growing both in axenic culture and 
within P. bursaria (Table 2). 

C. variabilis NC64A encodes at least 40 putatitve AA transporters 
(Table 1). Genes 1 to 35 in Table 1 were identified by Blanc et al. [4] 
on the basis of their significant homology with the AA transporter pro¬ 
tein family profile (Pfam PF01490.9). Of these 3 5 genes, 15 of them were 
expressed at significant levels as judged by read counts from RNAseq 


5 



NC64Agene ID 


Fig. 7. C. variabilis NC64A mRNAs coded by AA transporter genes during axenic growth. 
Normalized mRNA abundance of 15 AA transporter genes. 


experiments shown in Fig. 7 [4,42]. Each of these 35 genes had an appar¬ 
ent ortholog in C. sorokoniana UTEX-1230 (Table 1 ), and in 5 cases, two 
genes from NC64A identified the same gene from UTEX-1230 (Table 1 ). 

In a separate approach, the characterized AA transporters in the 
plant Arabidopsis thaliana (Table 1 in [46]) were used to perform recip¬ 
rocal BLAST searches against NC64A [4] and UTEX-1230 (UNL algal con¬ 
sortium, in preparation) genomes, using an expected value of 1 x 10“ 10 
as a cutoff. The major predicted isoform for each identified locus was se¬ 
lected and used to perform a BLAST search against the Arabidopsis ge¬ 
nome, again with an expected value of 1 x 10“ 10 as a cutoff. Each 
algal protein returned an Arabidopsis AA transporter, and the gene des¬ 
ignations and E-values for NC64A and UTEX1230 are presented in Sup¬ 
plementary Tables 1 and 2, respectively. Results of the initial and 
reciprocal BLAST searches identified 16 putative orthologs for 
Arabidopsis AA transporters in NC64A (Supplementary Table 1) and 11 
of these putative orthologs coincided with genes, which had been iden¬ 
tified by [4], see Table 1, while 5 of them had not been identified previ¬ 
ously. These 5 are indicated by genes 36 to 40 in Table 1. By contrast, 25 
putative orthologs were identified in UTEX 1230 (Supplementary 
Table 2). Due to the highly interrupted nature of the C. sorokoniana 
UTEX-1230 genome, there remain some discrepancies in the annota¬ 
tion. However, we are confident that even the partial nature of some 
UTEX-1230 AA transporter orthologs (Table 1 ) reflects AA transporters. 


Table 2 

Expression summary of predicted AA transporter genes in C. variabilis NC64A grown in culture or as a symbiont in Paramecium bursaria. Expression data for each gene represent manual 
counts of reads aligning to the corresponding genomic interval normalized per million mapped reads. 


NC64A Gene ID 

Protein ID 

Scaffold 

Axenic 

Symbiont 

Fold change 

Log 2 FC 

58128 

EFN54604 

13:78,798-84,676 (-) 

7 

99 

13.46 

3.7 

24724 

EFN54340 

14:418,634-420,275 (-) 

30 

8 

0.27 

-1.8 

36103 

EFN54400.1 

14:833,000-835,270 (-) 

3 

17 

6.01 

2.5 

53357 

EFN53996 

15:290,406-294,216 (-) 

30 

8 

0.27 

-1.8 

58448 

EFN53780.1 

16:626,917-630,844 (-) 

52 

6 

0.11 

-3 

138133 

EFN59501 

2:2,150,538-2,155,669 (-) 

16 

9 

0.57 

-0.8 

32765 

EFN51990 

25:222,544-226,977 (-) 

2635 

2794 

1.06 

0.08 

37093 

EFN51991 

25:227,152-230,758 (-) 

35 

42 

1.19 

0.26 

59057 

EFN51898 

25:231,079-235,117 ( + ) 



No Exp 


140447 

EFN58455.1 

3:553,777-556,173 ( + ) 

95 

32 

0.33 

-1.57 

142334 

EFN58316.1 

4:1,840,877-1,844,482 (-) 

48 

12 

0.25 

-1.95 

7483 

EFN50706.1 

43:40,210-42,680 ( + ) 

8 

8 

0.97 

-0.04 

17797 

EFN50713 

43:87,189-87,664 ( + ) 



No Exp 


57473 

EFN56726 

7:287,118-291,257 (-) 

83 

84 

1 

0 

59479 

EFN50622 

79:1541-9952 ( + ) 

1124 

4846 

4.31 

2.1 

144770 

EFN56324 

8:392,270-395,765 (-) 

25 

29 

1.14 

0.19 

AAT-A11 


Sum of all reads 

4193 

7996 

1.9 

0.93 
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3.9. Comparison of in symbiont vs. axenic expression 

The 16 AA transporters identified in NC64A (Supplementary Table 1) 
were examined for their expression levels in the axenic and endosymbi- 
ont states (Table 2). Fourteen of the AA transporter orthologs were 
expressed at detectable levels in the 2 P. bursaria assemblies (Table 2), 
and of these expressed isoforms, 3 were down regulated in axenic cul¬ 
ture, 6 were up regulated, and 5 remained roughly constant (Table 2). 
The sample size and differences in sequencing platform precluded a for¬ 
mal statistical analysis of the significance of these differential gene ex¬ 
pressions. However, 2 genes (32765 and 59479) accounted for 90 and 
95% of the mapped reads in axenic and symbiont growth, respectively; 
indicating that these transporters likely provide the majority of AA up¬ 
take capabilities in NC64A and other symbiotic Chlorella. The expression 
of 32765 is equivalently high in the axenic and symbiotic states while 
59479 is elevated 4.3-fold in the symbiont (Table 2). Thus, mining the 
genomes and transcriptomes available for Chlorella has led us to identify 
AA transporter genes that might play important roles in symbiotic N 
metabolism. It will be of interest to learn the AA specificity of the 
32765 and 59479 AAT genes. 

4. Discussion 

Nine Chlorella strains (5 symbiotic and 4 free-living) were compared 
with regard to their abilities to grow on organic and inorganic N sources. 
The 9 strains separated into two clusters, one containing all the symbi¬ 
otic strains and the other all the free-living strains. This analysis con¬ 
firms the robust nutritional/metabolic differences between symbiotic 
and free-living Chlorella strains as well as the generalization that the 
free-living Chlorella are better adapted for inorganic N sources while 
the symbiotic Chlorella are adapted for organic N sources. In particular, 
symbiotic Chlorella strains (a) could not use N0 3 ; (b) exhibited some 
NH 4 toxicity in that the cells grew slightly on 1 mM NH 4 but poorly or 
not at all on 10 mM NH 4 ; (c) in general preferred urea or amino acids; 
and (d) used two amino acids (Asn and Ser) that were not utilized by 
the free-living strains. We conclude that the symbiotic algae have phys¬ 
iological signatures that are conserved after they are separated from 
their symbiotic hosts and that one such signature concerns the constitu¬ 
tive expression of their AA transporters. For NC64A, it seems likely that 
these AA transporters reflect the mechanism by which P. bursaria might 
control the growth rate and population of its photosynthetic algal 
symbiont. 

Thus, we have physiological, genomic, and transcriptomic data 
showing that the symbiotic and free-living Chlorella strains differ signif¬ 
icantly in their N metabolism. The biological implications of this conclu¬ 
sion are three-fold. 

i) The inability of 4 symbiotic Chlorella to use N0 3 as the N source 
occurs despite the fact that recent work by Sanz-Luque et al. 
[43]) showed that C. variabilis NC64A has a complete set of 
genes needed for N0 3 assimilation (transporters, reductases, 
and synthesis of the molybdenum cofactor Moco). They conclud¬ 
ed that the mechanism responsible for silencing N0 3 utilization is 
unknown [43]. 


Our results provide a possible connection between the endosym- 
biotic life style, AA transporters, and virus susceptibility, illustrat¬ 
ing the trade-offs endosymbiotic Chlorella must make in nature, 
iii) Although the scope of this paper focuses mainly on N metabo¬ 
lism, the galactose preference exhibited by all symbiotic algae 
tested is also relevant. For instance, [4], suggested that NC64A 
has more enzymes involved in galactose metabolism than other 
sequenced chlorophytes. They looked for 2 protein families of ga¬ 
lactose related carbohydrate esterases, 18 families of glycosyl hy¬ 
drolases, and 8 families of glycosyl transferases present in the 6 
algal genomes which had been sequenced at that time. NC64A 
had a total of 149 galactose related genes whereas Chlamydomonas 
reinhardtii , Micromonas pusilla RCC299, Micromonas sp. CCMP, 
Ostreococcus lucimarinus, and 0. tauri had 84,36,30,36, and 29 re¬ 
spectively [4]. Thus genetic and physiological changes in the met¬ 
abolic capabilities of Chlorella symbiotic strains include major 
changes in their ability to assimilate both N and C sources. 

Future studies should improve our knowledge about the entirety 
of nutrient acquisition in Chlorella, including the assimilation of 
different N and C sources and its regulation, their relationship to 
the rest of metabolism, how metabolism is altered during symbio¬ 
sis, and how these aspects of symbiosis are tied to the viral suscep¬ 
tibility of some Chlorella strains. 


Supplementary data to this article can be found online at http://dx. 
doi.org/10.1016/j.algal.2016.06.001 . 
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